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We introduce a simple model of SO{N) spins with two-site interactions which is amenable to 
quantum Monte-Carlo studies without a sign problem on non-bipartite lattices. We present nu¬ 
merical results for this model on the two-dimensional triangular lattice where we find evidence for 
a spin nematic at small N, a valence-bond solid (VBS) at large N and a quantum spin liquid at 
intermediate N. By the introduction of a sign-free four-site interaction we uncover a rich phase 
diagram with evidence for both first-order and exotic continuous phase transitions. 


The destruction of magnetic order by quantum fluctu¬ 
ations in spin systems is frequently invoked as a route 
to exotic condensed matter physics such as spin liquid 
phases and novel quantum critical points [1-3] . The most 
commonly studied spin Hamiltonians have symmetries 
of the groups SO(3) and SU(2) which describe the rota¬ 
tional symmetry of 3-dimensional space. Motivated both 
by theoretical and experimental [4] interest, spin mod¬ 
els with larger-symmetries have been introduced, e.g. 
extensions of SU(2) to SU(A^) [5-8] or Sp(iV) [9]. 

The extension of SO(3) to SO(fV) is an independant 
large-fV enlargement of symmetry, with its own physical 
motivations [10]. While there have been many studies 
of SO(iV) spin models in one dimension [11-13], our un¬ 
derstanding of their ground states and quantum phase 
transitions in higher dimension is in its infancy. To this 
end, we introduce here a simple SO{N) spin model that 
surprisingly is sign free on any non-bipartite lattice. This 
model provides us with a new setting in which the de¬ 
struction of magnetic order can be studied in higher di¬ 
mensions using unbiased methods. As an example of in¬ 
terest, we present the results of a detailed study of the 
phase diagram of the our SO(A^) anti-ferromagnet on the 
two-dimensional triangular lattice. 

Models. - Consider a triangular lattice, each site of 
which has a Hilbert state of N states, we will denote the 
state of sitej as \a)j (1 < a < N). Define the A^(A^—l)/2 
generators of SO{N) on site i as with a < P; they 
will be chosen in the fundamental representation on all 
sites: = i6p.y\a)j — iSa-y\(3)j. Now consider the 

following SO{N) [14] symmetric lattice model for N > 3, 

= ( 1 ) 

(b> 

where the implies a summation over the N{N — l)/2 
generators and (ij) is the set of nearest neighbors. To see 
that Hj does not suffer from the sign problem, define a 
“singlet” state on a bond, l^y) = ^he 

singlet projector Using these operators 

and ignoring a constant shift we find the simple form [15], 

Hj = -JY.Pu- ( 2 ) 

(b> 


We make four observations: First, it is possible to cre¬ 
ate an SO(A^) spin singlet with only two spins for all 
N (in contrast to SU(A^) where N fundamental spins 
are required to create a singlet); Second Eq. (1) being 
a sum of projectors on this two-site singlet is the sim¬ 
plest SO(A^) coupling, despite it being a biquadratic in¬ 
teraction in the generators Third, since the singlet 
has a positive expansion, Hj is Marshall positive on any 
lattice; Fourth, on bipartite lattices Hj is equivalent to 
the familiar SU(N) anti-ferromagnet [6], i.e. the obvious 
SO{N) of Eq. (1) is enlarged to an SU(A^) symmetry. 
Since the bipartite SU (N) case has been studied in great 
detail in past work on various lattices [7, 16-23], we shall 
concern ourselves here with the non-bipartite SO(A^) case 
which is relatively unexplored. 

Phases of Hj: Starting at = 3, Eq. (1) becomes 

H = —^ J2{ij) with S the familiar S' = 1 rep¬ 

resentation of angular momentum. Previous numerical 
work has shown that this triangular lattice S = 1 bi¬ 
quadratic model [24, 25] has an SO(3) symmetry break¬ 
ing “spin nematic” magnetic ground state (we shall de¬ 
note this phase by SN). The ground state of Hj for N > 3 
has not been studied in the past. 

In the large-A^ limit, analogous to previous work for 
SU(A^) anti-ferromagnets on bipartite lattices [26], the 
ground state is infinitely degenerate and consists of dimer 
coverings where each dimer is in j Sij ). At leading order in 
1/N, Hj introduces off-diagonal moves which re-arrange 
parallel dimers around a plaquette, mapping Hj at large- 
A^ to a quantum dimer model on the triangular lattice 
with only a kinetic term, 

TJqdm = —t I f 1 —+ A.c.) (3) 

plaq 

where the sum on plaquettes includes all closed loops of 
length four on the triangular lattice. The ground state 
of this model has been found in previous analytic [27] 
and numerical work [28] to be a '/V2 x 1/12 valence bond 
solid (VBS), breaking the lattice translation symmetry. 
We thus expect that at large but finite values of N, Hj 
should restore its SO{N) symmetry and enter this same 
VBS state. 
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FIG. 1. Equal time structure factors for SN order [S'sN(k)], 
and susceptibility for VBS order [xvbs (k)] shown for = 10 
and N = 14, for the Hj model, Eq. (1) with L = 48. The 
Bragg peaks for SN (VBS) weaken (sharpen) with increasing 
N. The cartoon of the Brillouin zones shows the location of 
the ordering vectors of both order parameters. Quantitative 
finite size scaling of these orders is shown in Fig. 2. 

Since Hj has SN order for N = 3 and is expected to 
have a non-magnetic VBS at large-A^, it is interesting 
to ask what the nature of the transition at which SN 
magnetism is destroyed. The answer to this question is 
unclear based on current theoretical ideas and is best 
settled by unbiased numerical simulations. Exploiting 
that Hj has no sign problem we study it as a function 
of iV on L X L lattices at temperture /3 by unbiased 
stochastic series expansion [29] quantum Monte Carlo 
simulations, with a previously described algorithm 
[24]. The SN state is described by the matrix order 
parameter ~ The static structure 

factor, 5sN(k) = (Qaa(*)QaaO’)) 

is used to detect SN order. For the VBS order, we 
construct the k dependent susceptibility of dimer- 
dimer correlation functions in the usual way from 
imaginary time-displaced operators: xvBs(k) = 

^ / dr(Pr„r,+s(T)Pr,,r,+^(0)). 

Throughout this paper we have fixed j3 = L for our finite 
size scaling [15]. 

As shown in Fig. 1, a peak in S'sN(k) is found at the T 
point. Comparing the data at Af = 10 and N = 14, 
already qualitatively it is possible to see the peak in 
S'sN(k) softens as N is increased. In contrast xvBs(k) 
develops sharp peaks at the X and M points as N is 
increased. These are precisely the momenta at which 
previous numerical studies of the triangular lattice quan¬ 
tum dimer model Eq. (3) have observed Bragg peaks [28], 
validating the large-Af mapping to Eq. (3) made ear¬ 
lier. To detect at which N, the magnetic order is de¬ 
stroyed and the VBS order first sets in, we study the 
ratio, i?sN = 1 — (where a = x — y/V^) as 

a function of L. must diverge in a phase in which 

the Bragg peak height scales with volume and becomes 
infinitely sharp. On the other hand it must go to zero in 
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FIG. 2. Grossing plots of the ratios Rsn and i?vBS as a 
function of the discrete variable N for the Hj model. It is 
seen that spin nematic order is present for N < 10. VBS order 
on the other hand is present for A > 12. A = 12 appears 
to be on the verge of VBS order. Interestingly, A = 11 has 
no SN or VBS order. In the text, we present evidence that 
this phase is a QSL. The inset in the upper panel shows Rsn 
scales to zero at A = 11, despite non-monotonic behavior at 
intermediate L. 

a phase in which the correlation length is finite and the 
height and width of the Bragg peak saturate with sys¬ 
tem size. At a critical point standard finite size scaling 
arguments imply that the ratio, i?sN becomes volume in¬ 
dependent. All of these facts together imply a crossing 
in this quantity for different L. Fig. 2 shows the i?sN 
and i?vBS ratios (an analogous quantity constructed for 
the VBS order from xvBs(k) close to the M-point) as a 
function of the discrete variable N for different L. The 
data for i?sN shows that the magnetic order is present 
for N < 10. The i?vBS data shows that the long-range 
VBS order is present for N > 12. From Fig. 2 we find 
that A = 12 is on the verge of developing VBS order; 
from the system sizes accesible we are unable to reliably 
conclude whether A = 12 has long range VBS order or 
not from our study. However, taken together the data 
show definitively that A = 11 has neither VBS nor SN 
order. As we shall substantiate below, at A = 11, Hj is 
a quantum spin-liquid (QSL). 

J-Q models: In order to clarify the global phase dia¬ 
gram of SO (A) anti-ferromagnets and access the quan¬ 
tum phase transitions between the SN, VBS and QSL 
phases found in Aj, it is of interest to find an interaction 
that can tune between these phases at fixed A. In order 
to be meaningful, the new coupling must preserve all the 
symmetries of Hj. To this end, we introduce and study a 
generalization of the four-site Q term of SU(2) spins [30], 

Hq = -QY1 (4-Pfe; + PaPjk) (4) 

{ijkl) 

where the sum includes elementary plaquettes of length 
four on the triangular lattice (with periodic boundary 
















3 




FIG. 3. Phase diagram of jfljQ [Eqs. (2,4)] for different 
values of N. The left panel shows the phase diagram for 
small N, where a first order SN-VBS transition is found for 
G < N < 9, (see Fig. 4). As N is increased we find the 
first order transition weakens. The right panel shows how an 
intermediate QSL phase emerges for = 10 and = 11. 
Transitions from the QSL to both SN and VBS phases are 
continuous on the large systems studied, see Fig. 5. 

conditions on an L x L system there are 3L^ such pla- 
quettes). For a fixed-Af, Hq provides a tuning parameter 
which preserve both the internal and lattice symmetries 
of Hj and hence allows us to study the generic phase 
diagram of SO(Af) magnets. A summary of the phase 
diagram of Hjq in the N-Q/ J plane is in Fig. 3: The Q- 
interaction destroys the SN order and gives way to VBS 
order only for N > 6. We have found evidence for direct 
first-order SN-VBS transitions for 6 < V < 10 and exotic 
continuous SN-VBS transitions for N = 10 and N = 11. 

As an example of our observed first-order behavior we 
present in Fig. 4, our study of the N = 7 QMC data for 
the spin stiffness ps = {W^) jL (where Wx is the winding 
number of the spin world lines), which acts as a sensitive 
order parameter for the SN phase, and the VBS order 
parameter Oyes — XVBs(M)/Afsite- Clear evidence for a 
direct first order SN-VBS transition at V = 7 is found. 

The nature of the transition changes at N = 10, where 
evidence for two phase transitions is found. As shown 
in Fig. 5 the SN order vanishes at a Q/J smaller than 
the value at which VBS order develops. Although the 
difference is small for N = 10, it is significant. The data 
for V = 11 in Fig. 5 shows that the SN and VBS orders do 
not vanish at the same point. In fact i?sN indicates that 
the SN order has vanished already at Q/J = 0, consistent 
with our previous analysis of IIj. As illustrated by the 
dashed and solid lines in Fig. 3, the appearance of the 
QSL phase is consistent with a global phase diagram for 
the SO(N) magnets. 

QSL phase and criticality: We have identified the 
ground state between SN and VBS as a QSL, since it 
does not show evidence for any Landau-order. Were the 
intermediate phase characterized by a conventional or¬ 
der parameter, we would have expected strong first order 
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FIG. 4. First-order SN to VBS transition in Hjq at N = 7. 
The upper panel shows the VBS order parameter and the 
stiffness as a function for Q/J for different L indicating a di¬ 
rect SN-VBS transition. The lower panel shows MC histories 
(and histograms in the inset) at Q/J = 1.26, providing clear 
evidence that the SN-VBS transition at A = 7 is direct and 
first-order. 

transitions of the kind between SN and VBS (see Fig. 4), 
instead we find continuous transitions. 

There are field theoretic reasons to expect a QSL on 
quantum disordering a spin nematic. The long-distance 
description of our SO(V) models is given by a 
theory (in contrast to the description of SU(V) 

models [31]), which can be described as N real matter 
fields coupled to a Z 2 gauge field. Such a theory is ex¬ 
pected to host three phases [32], a symmetry breaking 
phase in which the matter condenses (which we identify 
in our spin model as the SN), a stable phase in which the 
matter gets a gap and the Z 2 gauge theory is deconfined 
(identified here as the QSL) and a phase in which matter 
is gapped and the Z 2 is confined (identified here as the 
VBS). Thus, the SN-QSL critical point should be in the 
universality class of 0(A)* critical point [3]. The QSL- 
VBS phase transition should be in the same universality 
class as the critical point between these identical phases 
in the quantum dimer model since the magnetic fluctu¬ 
ations are gapped in both the QSL and VBS phases. A 
previous analysis of this phase transition has predicted 
an 0(4)* phase transition [27], where the VBS order pa¬ 
rameter is identified with a bilinear of the primary field. 

A detailed study of the critical phenomena at A = 10 
and A = 11 is clearly beyond the scope of the current 
manuscript. We shall be satisfied here with a brief anal¬ 
ysis: At the QSL-VBS critical point, we are able to carry 
out reasonable data collapses [15] at both A = 10 and 
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FIG. 5. Crossings of Rsn (above) and Rvbs (below) sig¬ 
naling the location of the onset of long-range SN and VBS 
orders at A'" = 10 (left) and = 11 (right). At A'" = 10, 
7?sn and Rvbs cross at close but significantly different cou¬ 
plings, Qc = 0.100(5) and Qc = 0.117(2) respectively. At 
A = 11, Rsn appears to have crossed &t QjJ <0 (we can¬ 
not study this region because of the sign problem), whereas 
Rvbs crosses at Qc = 0.042(3). From the location of the 
crossings, for both A = 10 and A = 11, we can infer an in¬ 
termediate phase which is neither SN nor VBS, as shown in 
Fig. 3(b). We present arguments that this phase is a QSL. No 
direct evidence for first order behavior is found at either of 
the transitions, though a weakly first order SN-QSL cannot 
be ruled out. The QSL-VBS transitions shows good scaling 
behavior with unconventional critical exponents. 


N = 11 for OyBS both X and M ordering vectors, 
see Fig. 1) and i?vBS) where we find, ryvBS = 1-3(2) and 
iAbs = 0.65(20) for the anomalous dimension of Ovbs- 
The unusually large value of r?vBS is a direct consequence 
of fractionalization in the intermediate QSL phase and 
is often regarded as a smoking gun diagnostic of ex¬ 
otic critical points (see e.g., [33]). More quantitatively, 
our critical exponents are in rough agreement with the 
best estimate of r; = 1.375(5) of the bilinear field and 
ly = 0.7525(10) in the 0(4) model [34]. We note that the 
values for jyvBS and agree within the quoted errors 
for A = 10 and A = 11. Taken together, this bolsters 
the case that the intermediate QSL phase has Z 2 frac¬ 
tionalization, albeit more work is needed for a definitive 
identification. Unfortunately, the SN-QSL transition, ob¬ 
served only at A = 10, has large corrections to scaling 
and we are unable to reliably determine its critical ex¬ 
ponents or determine whether it is a weakly first order 
transition (no direct evidence for a first-order transition 
has been found of the type shown for the A = 7 case). 

In summary, we have introduced a new family of sign- 
free SO (A) spin models, which can be regarded as non- 
bipartite generalizations of their popular SU(A) cousins. 
The triangular lattice model which we have studied thor¬ 


oughly here hosts a spin nematic, a VBS with a large unit 
cell, a quantum spin liquid phase and unusual quantum 
critical points. The absence in the SO(A) models of a di¬ 
rect continuous “deconfined quantum critical point” [33] 
is in striking contrast to previous simulations of the re¬ 
lated bipartite SU(A) models [8, 23]. We have offered a 
plausible field theoretic scenario that naturally explains 
this difference. It is interesting that the absence (pres¬ 
ence) of a QSL in bipartite SU(A) (non-bipartite SO(A)) 
spin models seems to track the absence or presence of 
this phase in the kind of quantum dimer models that our 
model maps to at large-A [35]. 

While the study in this paper has focussed on the tri¬ 
angular lattice, our family of models, Eq. (2,4) may be 
constructed sign free on any two or three dimensional 
non-bipartite lattice. Because of the larger degree of frus¬ 
tration, the kagome system may provide a wider swath of 
the QSL phase and hence could possibly allow a more de¬ 
tailed study of this phase, even if the phase diagram is of 
the same form found here. Exploring the phase diagram 
and quantum phase transitions of the three dimensional 
pyrochlore system is an exciting open direction for future 
work. 
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SUPPLEMENTARY MATERIALS 
Model and Symmetries 

Here we provide some additional details of the models 
introduced in Eq. (1) and Eq. (2). 

Mapping between Eqs. (1) and (2) 

To see the connection between the two Hamiltoni¬ 
ans Eq. (1) and Eq. (2). We consider two SO{N) 
spins. We can combine them into three representa¬ 
tions: a singlet (S), symmetric (y) and anti-symmetric 
(<&) representations of dimensions: 1, -I- ^ — 1 and 

^ - Y- Now construct projectors on these representa¬ 
tions, Ps,Px Clearly Ps + P^ + P^ = 1 and 

Pg = PsjP^ = -PjcjC’l = It is straightforward to 
show that, {Li ■ Lj) = —{N — l)Ps — P^ + P^ by explic¬ 
itly acting on the symmetrized wave-functions. From this 
it follows that {Li ■ Lj)^ = (TV — l)^Ps + P^P P<s>- From 
which it follows that {Li ■ Lj)^ = ((TV — 1)^ — l) Pg -|- 1, 
which proves as claimed that for TV > 3, Eq. (1) and 
Eq. (2) are equivalent up to a constant. 

N = 2 

Although not studied in this manuscript, for the sake 
of completeness, we discuss our model at TV = 2. Even 
though Eq. (1) is trivial for TV = 2 (since squaring 
the only SO(2) generator is just an identity operator), 
Eq. (2) is a well defined non-trivial model. Identify¬ 
ing the two colors with f and | spins, Eq. (2) becomes 
H = E(^J) -JiSfSJ + SfSf) + JSfSl Previous work 
on the triangular lattice TV = 2 model has found clear 
evidence for SO(2) symmetry breaking superfluid order 
[1-3]. 

Symmetries 

We now discuss the symmetries of the model Eq. (2). 
We begin by observing that this model is invariant under 
uniform 0(TV) rotations where we multiply each basis 
state by an orthogonal matrix (one that satisfies 0^0 = 
I), since this leaves the singlet state invariant, i.e. 

Oa>yOa.7]\'^Tl') — ^|aa). (5) 

O' a a 

However we should identify rotations that only differ by 
changing all the local basis states by the same phase (in 
this case a sign). Here it becomes necessary to distinguish 
between even and odd TV. This is because the matrix — 1 
has determinant 1 for even TV and -1 for odd TV. Thus 


size 

N 

Q 

/3qmc 

Eex 

Eqmc 

2x2 

4 

0 

16 

-1.5 

-1.49997(3) 

2x2 

4 

1 

16 

-4.5 

-4.5000(1) 

2x2 

5 

0 

16 

-1.4 

-1.40000(5) 

2x2 

5 

1 

16 

-4.2 

-4.2001(1) 

2x3 

3 

2 

16 

-6.0657499233 

-6.0656(1) 

3x3 

2 

0 

16 

-1.7026987262 

-1.70269(2) 

3x3 

2 

1 

16 

-3.7290340614 

-3.72900(3) 

3x3 

2 

2 

16 

-5.7639923092 

-5.76402(5) 


TABLE I. Test comparisons of ground state energies from ex¬ 
act diagonalization and average energies from finite-T QMC 
studies of the SO{N) model introduced here. Note that J = 1 
always. The energies reported here are per site and on tri¬ 
angular lattices with periodic boundary conditions such that 
there are always 3L^ bonds and 3L^ plaquettes in Eqs. (2,4). 
This causes some terms to be appear more than once for the 
2x2 and 2x3 systems. 

for odd-TV the symmetry is simply SO (TV), since the rest 
of 0(TV) is obtained from SO(TV) by multiplying by -I. 
For even-TV however SO (TV) has pairs of elements that 
cause the same basis transformation up to a sign, e.g. 
1 and -1. On the other hand unlike the case of odd-TV, 
the 0(TV) matrices with determinant -1 are independant 

symmetries, so the symmetry realized for even-TV is an 
0(N) 

Z2 ■ 

Symmetry on Bipartite Lattices 

On bipartite lattices the orthogonal rotation symme¬ 
try, Eq. (5) gets extended to a unitary symmetry (with 
U^U = I) so long as the singlet is defined between sites 
on opposite sub-lattices, and A sub-lattice spins are ro¬ 
tated by U and B sub-lattice spins are rotated by U*, 

^\aa) ^^U*^Uar,hri) =^\aa). ( 6 ) 

a Q. oc 

Since a uniform phase change of the all the states locally 
does not have physical consequences, the model is said 
to have an SU(TV) symmetry, as has been discussed and 
extensively studied previously in such models, see e.g. 
Ref. [4] for a review. 

Ground state theorems 

Marshall’s sign theorem guarantees that the ground 
state of Hjq = Hj + Hq is an SO(TV) singlet. In addi¬ 
tion, on the triangular lattice, which is the focus of our 
study here, there is no simple translationally invariant 
covering of two-site singlets, leading us to suspect that 
a generalization of the SU(2) square lattice Lieb-Schultz 
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FIG. 6. Estimation of optimal /3/L ratio for our simulations 
by comparison of the fluctuations of the temporal ((W^)) and 
spatial {(Wx)) winding numbers . Data shown is for Hj at 
N = 10. 


1.42) and that the crossing point does not move much 
with system size. Thus for simplicity we have chosen 
P/L = 1 throughout the paper. 




Mattis (LSM) theorem [5] applies to Hjq on this lat¬ 
tice, i.e. in the thermodynamic limit there must be a 
degeneracy in the ground state, so that a simple gapped 
paramagnet is not possible - either a symmetry is bro¬ 
ken or the ground state is exotic. A rigorous proof of 
this intuitive assertion is expected to be at least as tech¬ 
nical as the proof for the bipartite N = 2 case [5] and 
is beyond the scope of this work. As we saw above, on 
a one-dimensional chain, which is bipartite, our model is 
equivalent to the SU(A^) model studied by Affleck [6] and 
is hence expected to have an LSM degeneracy. 

Numerical Simulations 

QMC energy tests 

Here we provide the results of some QMC tests on small 
lattices for the total energy per spin of our models, Hj + 
Hq, Eqs. (2,4) on the triangular lattice, for completeness 
and future comparisons. 

Choice of P = L 

In Fig. 6 we show the dependence of the fluctuations of 
the temporal and spatial winding numbers on the “aspect 
ratio” of our simulation cell. /3 is a measure of the extent 
of the imaginary time and L is an estimate for the linear 
spatial extent. We study the fluctuations of the temporal 
and spatial winding numbers as the ratio j3/L is varied 
for two different sizes, L = 32 and L = 48 at = 10 in 
the model Hj. We find that both quantities are balanced 
at a value of /3/L which is of the order of one (close to 


FIG. 7. Finite size scaling of spin nematic order parameter 
and spin stiffness for the model Hj, in Eq. (1) with 10 < A < 
14. Shown on the left is the square of the spin nematic order 
parameter, 0 |n (the height of the Bragg peak in Fig. 1), and 
on the right is the spin stiffness, ps, plotted as a function 
of 1/L. The data confirms that the system is magnetically 
ordered for A < 10 and non-magnetic for A > 11. 



FIG. 8. Finite size scaling of VBS order parameter for 
the model Hj, in Eq. (1) with 10 < A < 14. For A > 13 
we encounter difficulties in equilibrating the system for sizes 
larger than L = 48 due to formation of long-range VBS order. 
Rather than extrapolating Oyssi study the crossing of the 
ratio i?vBS, which provides a reliable way to detect the onset 
of long-range VBS order on moderate system sizes [see Fig. 2]. 














Extrapolation of order parameters 


-Rvbs.m 


The simplest estimate for long range order is to study 
whether the height of the Bragg peak per unit volume 
extrapolates to a finite quantity in the thermodynamic 
limit. Unfortunately, this method becomes increasingly 
unreliable when the measured order is weak, e.g., close 
to a critical point. In such cases, results from extrapola¬ 
tions will depend on the form of the extrapolation used. 
A thorough discussion of these difficulties in quantum 
spin systems may be found in the literature [7]. It is 
for this reason that we prefer to work with the R ra¬ 
tios defined in the text. The disadvantage is that we 
do not know the order paramater in the thermodynamic 
limit, but the advantage is we can calculate the critical 
coupling reliably by studying the crossing of the R ra¬ 
tio. For completeness we present here the data required 
for extrapolation of both SN and VBS order parameters 
for Hj. To test quantitatively for long range order we 
study the scaling of the height of the peak in S'sN(k;), 
Ogj^ = S'M(k; = 0)/Aljite and the spin stiffness ps on fi¬ 
nite size systems with A^site = Lx L. Both quantities are 
expected to be finite in the M state and zero when the 
0{N) symmetry is restored. Fig. 7 shows finite size data 
for both quantities for different values of N. From these 
plots we conclude that the M symmetry is broken up to 
= 10 and is restored for A^ > 11, because scales to 
zero for these N. This behavior is mirrored in ps, albeit 
for intermediate L there is some non-monotonic behavior 
for A^ = 11. This is consistent with our conclusions in 
the main text made from the analysis of i?sN 

Finite size scaling for the VBS order parameter is 
shown in Fig. 8. Notice for the cases where there is VBS 
order {N = 13,14) we only have data for L < 48. For 
system sizes larger than this we face serious equilibra¬ 
tion issues with QMC as is expected, since the simula¬ 
tion gets locked into a symmetry broken VBS state. The 
plot serves to illustrate te ambiguity faced by making di¬ 
rect extrapolations. On the other hand, a study of the R 
ratios shown in Fig. 2 provides a more clear cut way to 
locate the critical point. 


QSL-VBS and SN-QSL phase transitions 



2.0 


1.5 


1.0 


0.5 


-20 0 20 
(Q - OJiV- 



FIG. 9. Examples of our data collapses at QSL-VBS critical 
point at A = 10 shown here with the best parameters deter¬ 
mined individually for each observable. The parameters used 
are Qc = 0.1170,0.1187,0.1171, 1/iwbs = 2.171,2.152,1.733 
and ?7vBS = 1.414,1.106, —. 



(Q - Q,)LL- {Q - Q^)LL-' {Q - q^)l^l 

FIG. 10. Examples of our data collapses at QSL-VBS 
critical point at A = 11 shown here with the best param¬ 
eters determined individually for each observable. The pa¬ 
rameters used are Qc = 0.04330,0.04182,0.03917, 1/vybs = 
1.902,1.858,1.244 and t^vbs = 1.383,1.133, -. 


Here we present some details of the study of the both 
the QSL-VBS (A = 10,11) and SN-QSL (A = 10) phase 
transitions found in our model. 

We obtain critical exponents at the QSL-VBS critical 
point by attempting a data collapse, see Fig. 9,10. We 
use the standard finite size scaling ansatz for the order 
parameter and the crossing ratio, 

(OvBs) = (7) 

= ( 8 ) 


where g = {Q — Qc)/ J- We continue to work with (3 = L 
as discussed. No attempt is made to make use of correc¬ 
tions to this leading scaling behavior. Our main objec¬ 
tive is to determine the universal number ryvBS for the 
QSL-VBS transition for A = 10 and A = 11. We find 
acceptable collapses for our data sets over a wide range 
of ^Vbs- On the other hand, the estimate for tjybs is 
relatively stable over our various fits. The values and er¬ 
rors of the critical exponents quoted in the main text are 
based on the variation observed by using different data 
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FIG. 11. Crossing of L and Lj^ for various dimensionless 
quantities for N = 10. The dashed line is a quantity (Rvbs) 
whose crossing locates the VBS transitions and the solid lines 
are for two independent quantities {ps and Rsn) that locate 
the SN transition. The blue and red semi-circle shows the 
range of critical couplings based on the extrapolation of this 
data. Note that (1) there is clear evidence for an intermediate 
phase. (2) the errors in the SN transition are signihcantly 
larger than those for the VBS transition. 


sets. A higher precision study should be possible with 
access to more accurate data and larger system sizes. In 
order to carry out the collapse numerically, we make use 
of a recently developed Bayesian approach to scaling [8] . 
We note that difficulties in obtaining accurate values of 
the critical exponents at exotic transitions in quantum 
spin models is a well-documented difhculty [9] . 

Another quantum phase transition takes place between 
SN and QSL. In our model this transition appears only at 
N = 10. In Fig. II we study the drift of various crossing 


quantities at the critical points. Presumably the signifi¬ 
cant drift for the crossing at the SN-QSL transition are 
due to corrections to scaling. We have looked for signs 
of first order behavior as we found for smaller-and not 
found them here, though the possibility of a very weak 
first order transition cannot be ruled out. The corrections 
to scaling hamper efforts to extract critical exponents at 
this phase transition. In contrast the QSL-VBS transi¬ 
tion shows a reasonably converged crossing point with 
a nice scaling regime, where the crossing points do not 
depend significantly on L. 
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